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Abstract. Using the exact computation of a large number of moments of the 
distribution function of the end-to-end distance G{r^ N) of the worm-like chain, we 
have established the analytical form of the coefficients in Taylor expansions of the 
moments for short chain lengths N . The knowledge of these coefficients enabled us 
to resummate the moment expansion of G'(r, N) by taking into account consecutively 
the deviations of the moments from their stiff rod limit. Within this procedure we 
have derived the short-chain expansion for G{r,N), the scattering function, and the 
extension-force relation, which take into account the deviations of the moments from 
their stiff rod limit to the seventh order in N. 
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1. Introduction 

Polymers with contour lengtli L mucli larger than the persistence length /p, which is the 
correlation length for the tangent-tangent correlation function along the polymer and 
is quantitative measure of the polymer stiffness, are flexible and are described by using 
the tools of quantum mechanics and quantum field theory PJ-IH]- If the chain length 
decreases, the chain stiffness becomes an important factor. Many polymer molecules 
have internal stiffness and cannot be modelled by the model of flexible polymers 
developed by Edwards p. 

The standard coarse-graining model of a worm-like polymer was proposed by 
Kratky and Porod [7j. The essential ingredients of this model is the penalty for bending 
energy, and the local inextensibility. The latter make the treatment of the model much 
more difficult. There is a substantial amount of studies of the Kratky-Porod model 
in the last half century [H-JHl (and the references therein). In recent years there is 
an increasing interest in theoretical description of semifiexible polymers jTHj-ISni (and 
the references therein). A reason for this interest is due to potential applications of 
semifiexible polymers in biology and in researches on semi- crystalline polymers. 

In this article we present results of the study of the behaviour of the worm-like chain 
for short lengths. The consideration is based on our recent work |21],[^, where the 
Fourier-Laplace transform G{k,p) of the end-to-end distribution function G{r,N) {p is 
the Laplace conjugate to A^) was represented as the matrix element of the infinite order 
matrix P{k,p) = {I+ikDM)^^D with matrices D and M related to the spectrum of the 
quantum rigid rotator. A truncation of P{k, p) by matrix of order n gives the end-to-end 
distribution function G{k,p) as a rational function being an infinite series in powers of 
k'^, i.e. it contains all moments of the end-to-end distribution function, and describes 
the first 2n — 2 moments exactly. In context of eigenstates of the quantum rigid rotator, 
the truncation at order n takes into account the eigenstates with quantum number of 
the angular momentum up to the value I = n — 1. The moment expansion of G{k,p) 
can be represented as double series in powers of k"^ and 1/p. Alternatively, one can 
expand G{k, N) in double series in powers of k"^ and A^. Analyzing the series of known 
moments (we have analytically calculated the first 50 moments), we have established 
that the coefficients in subseries in powers of fc^ have a simple analytical structure, which 
enables one to perform a partial resummation of the moment expansion of G{k,p) or 
G{k,N) (see below). This resummation procedure results in an expansion of G{r,N) 
around the stiff rod limit, where the end-to-end distribution function is given by the 
expression, GT:{r,N) = l/{4'n'N'^)6{N — r), so that because G'r(r, A^) is a distribution, 
so derived short-chain expansion is not a Taylor expansion but rather an expansion in 
the space of distributions. The knowledge of G{k, N) enables one to compute directly 
the scattering function S{k, N) and the extension-force relation R{f)- In this paper we 
present the results of the calculation of terms of the short-chain expansions of G{r, N), 
S{k, N), and R{f) by taking into account the deviations of the moments from their stiff 
rod behaviour to the seven order in chain length A^. The procedure can be extended to 
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higher orders. The short-chain expansions besides the interest in their own can be used 
for comparisons with approximative treatments, and also in studies of the behaviour of 
short semiflexible polymers. 

The present paper is organized as follows. Section |21 introduces to the description 
of the worm-like chain using the formalism of quantum rigid rotator. Section El explains 
the idea of the derivation of the short-chain expansion, and presents the short-chain 
expansions for the distribution function of the end-to-end distance, scattering function, 
and extension-force relation. 



2. The formalism 



The Fourier transform of the distribution function of the end-to-end polymer distance 
of the continuous Kratky-Porod model [Zj G(k, L) = / (i^i?exp(— ik(R — Ro))G'(R — 
Ro, L) is expressed by the path integral as follows 

G(k,L) = j Dtis)^5{tisf-l)eM-^^jystis)-^j£ds{^Y), (1) 

where Ip is the persistence length, and t(s) = dr{s)/ds is the tangent vector at the point 
s along the contour of the polymer. The product over s in Eq.((T)) takes into account that 
the polymer chain is locally inextensible. In the following the arc length of the polymer 
L will be measured in units of Ip and will be denoted by A^. We now will consider the 
Green's function P{9,ip,N]9o,ipo,0) associated with Eq.dTJ. The differential equation 
for P is 

-^Pie, ^, N; ^0, ^0, 0) - l^l^p + u{n)p = siN)sin - n^), (2) 

where U(ktfi) = zktfj is the potential energy of the rigid rotator in an external field ik, 
where k is measured in units of The end-to-end distribution function is obtained 
from P{Q, N; Qo, 0) as follows 

G{k,N) = ^Jdnj dnoP{n,N;no,o). (3) 

The differential equation (jSJ can be rewritten as an integral equation as follows 

P{n, N; Qo, 0) = Po{n, N; Qq, 0) - ^ ds j dfl'Po{n, N; Q', s)U {ktn') P {Q' , s; Qq, 0), (4) 

where the bare Green's function Po{9, Lp, N; 9q, ipo, 0) reads 

Po(g, V, N; 00, v^o, 0) = ^ exp(- ^^^ ^^^^^ )YUO, ^1%, ^o), (5) 

with Yim{0,if) being the spherical harmonics, and / and m are the quantum numbers 
of the angular momentum. Due to the convolution character of the expression (0]) 
with respect to the integration over the contour length {Po{fl, N; Q', s) depends on the 
difference N — s), the Laplace transform of A^; flo, 0) in Eq.(@]) with respect to N 
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permits to get rid of integrations over the contour length. Thus, in the following we will 
consider the Laplace transform of G{k, N) with respect to A^. 

It was shown in |21],|2nj that the solution of Eq.(^ results in the following 
expression of the Fourier-Laplace transform of the end-to-end distribution function as 
the matrix element of an infinite order square matrix 



G{k,p) = {0\P'{k,p)\0) (6) 

with 

P'{k,p) = (/ + ikDM'y^D, (7) 
where the square matrices and D are defined by 

Mil, = + wi+i6i+i^i> , (8) 



with wi = ^P/{AP - 1), and 

respectively. The superscript s specifies that the quantities P'^ and ikf are square 
matrices. The quantity (0 | P^{k,p) \ 0) denotes the (1, 1) matrix element of the infinite 
order square matrix P'^. Since the (1, 1) matrix element corresponds to the expectation 
value of the quantum rigid rotator in the ground state with the quantum number / = 0, 
we prefer to use the above notation. Summations over the magnetic quantum number 
in the intermediate states in the expression of G{k,p) obtained by solving Eq.(jll) can 
be eliminated [201, so that summations in the intermediate states occur only over the 
eigenvalues of the angular momentum / = 0, 1, .... This is the reason why the calculation 
of G{k,p) reduces to the computation of the matrix element of an infinite order square 
matrix. 

The quantity P'^{k,p) plays the key role in the theory similar to the bare propagator 
in common quantum field theories. The end-to-end distribution function G{k,p) is 
simply the matrix element (^0\P^{k,p)\0^, the scattering function of the polymer is the 
inverse Laplace transform of G{k,p)/p'^ multiplied by 2/N, the partition function of the 
stretched polymer is Z{f,N) = G{k = -if/kBT,N) etc gHl- 

3. The behaviour of short worm-like polymer 
3.1. The idea of the short-chain expansion 

To derive the short-chain expansion for the Kratky-Porod chain we insert the expansion 
of D given by Eq.Q according to 

1 _ 1 (-1)" (/(/ + !))" 1 
^Kl + l)+P ~P^o ^ 
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into the moment expansion of G{k,p) 

-. oo 

G{k,p) = - J^^l-l)" (0 I {DM'f^ I 0) (A;^)'", (10) 

^ m=0 

and obtain as a result G{k,p) as double series in powers of {k/pY and 1/p as follows 

oo cx> /;2\" 



s=l ^ n=0 



The analysis of 50 exactly computed moments of G{r,p), which are proportional to the 
coefficients (0 | {DM'^y"^ \ 0) in Eq. lfTmi . shows that the coefficients are polynomials 
in powers of n of the order 2s— 3, i.e. they have the form = alra+a2^^^ + ...+a2s-3^^-^*~^ 
(s > 2). We determine the coefficients at given s using 2s — 3 known terms of 
the moment expansion. Further we test the correctness of so obtained from the 
comparison with the remaining known terms in the moment expansion, which were not 
used to determine Oj. The exact knowledge of 50 moments enables one to determine 
for s < 11. The coefficients we have obtained in this way are 

1 



^1 

c: 



n 



2n + 1' 

n 

~3' 
1 



Cl = -n{n + l){lAn + l), 

Cl = -^—n{n + l){n + 2){2n + 3)(508n3 - SAn"^ - 19n + 15), 
37800 

^" = - 3742200 ''^'' + 1)(^ + 2)(2n + 3)(2n + 5)(10220n^ - 6236n3 
+ 1597^2 + 737/2 - 372), 



We have also computed the coefficients Cj^ and C^, but do not write them here to 
save space. The determination of higher coefficients demands the knowledge of more 
moments of the end-to-end distribution function and can be obtained in a similar way. 
After the determination of the coefficients we have checked the equivalence of ()lip 
with (fTUI) to the corresponding order. Unfortunately, we did not succeed in deriving the 
expression for the coefficients for arbitrary s. 

With known coefficients one can sum the series over n in (fTT|). Restricting 
ourselves to the leading order Smax = 1 (smax is the number of terms in the sum over s 
in dni)) we obtain G{k,N) 

and the distribution function of the end-to-end distance 
a(r,iV) = ^5(iV-r). 
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in the stiff rod limit. Talcing into account the next term in the sum over s in Eq.flllj) 
gives the correction to the stiff rod limit of the end-to-end distribution function. Thus, 
the resummation of the moment expansion of G{k,p) according to (jllj) yields the short- 
chain expansion for Kratky-Porod chain. 

The above derivation does not allow to make claims on the convergence of the 
short-chain expansion. The comparison of results of computations of quantities under 
consideration (for example the scattering function or deformation-force relation) for 
different Smax gives a criterion determining the quality of the short-chain expansion. 

Notice that alternatively one can first carry out the inverse Laplace transform of 
and then expand the moments in Taylor series in powers of A^. As a result one will 
arrive at double series similar to (fTT|) with 1/p replaced by A^. This makes more clear 
the meaning of the resummation procedure. To leading order (smax = 1) one replaces 
the moments by its stiff rod behaviour, (i?^") = A^^". In the next to leading order 
(smax = 2) one takes into account the next-order corrections to the stiff rod behaviour 
of all moments, and so on. 

Due to the fact that the scattering function of the semifiexible polymer S{k, N) is 
the inverse Laplace transform of G{k,p)/p'^ multiplied with the factor 2/N, the short- 
chain expansion of G{k,p) enables one to get the short-chain expansion of the scattering 
function. In fact, Hermans and UUman |H| derived the stiff rod limit of the scattering 
function using the stiff rod limit of the moments. 

In following subsections we will consider separately the short-chain expansion of 
the end-to-end distribution function, the scattering function, and the extension-force 
relation, which also can be obtained from G{k, N). 



3.2. Distribution function of the end-to-end distance 

The series over n in (fTTjl for known coefficients can be easily expressed through the 
derivatives of the geometric series. Carrying out the inverse Laplace transform over p 
gives the short-chain expansion of G{k, N). We have computed G{k, N) by taking into 
account eight terms in the sum over s in (fTT|) i.e. Smax = 8. To save place we present 
below the result for Smax = 5 

G,{k, N) = G,{k, AT) + i ^^^^^ - cos (kN) 

+— ^^^(^^) +J_ksm (kN) N^- — cos (kN) + cos (kN) N'k^ 

630 k 5040 ^ ^ 630 ^ ^ 15120 ^ ^ 

1 Bin (*iV)iV3 1 m i3,i„(,^)^r 



5040 k 5600 ' ' 604800 

1 ^3 

cos (kN) cos (kN) (12) 



5040 ' ' 151200 
The subscript at G and at quantities below is Smax- Since Gs{k,N) tends to one for 
k 0, the end-to-end distribution function is normalized. It is easily to see that the 
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terms in (jl2|) can be represented as derivatives of sin {kN) jk with respect to i.e. of 
NGj-lk, N). The short-chain expansion of the distribution function of the end-to-end 
distance is then immediately obtained from ()12|1 as 

nG.ir, N) = ±6 {N - r) + ^6 {N - r) - l^^^) {N - r) + {N - r) 

r 7 9 31 

(N-r) — r25(2) m ^r) + —r6(N -r) - -—r^ 5^^'^ (N - r) 

30 ^ ' 1440 ^ ^ 63 ^ ' 1008 ^ ' 

■ rH^^^ {N-r)- r^5(3) (AT - r) + 47 {N-r)- -^r5^^^ {N - r) 



1440 ' ' 60480 ' ' 144 ' ' 1008 

where ^^'^•'(x) denotes the kth derivative of the Dirac's deha-function. Note that G{r, N) 
is a distribution function with respect to r, while the contour length is a parameter. 
This should be taken into account in using Eq. (|T!?|) to compute different mean values, for 
example the moments of the end-to-end distribution function of a short chain. Similar to 
the corresponding expansion of the function Sa{x) = l/(27ra^)^/^ exp(— a;^/2a^) in powers 
of a, the short-chain expansion of G{r, N) cannot be interpreted as a Taylor series but is 
rather a generalized Taylor expansion in space of distributions. The expansions p2H13p 
can be used to test approximative expressions of the end-to-end distribution function. 
We did not succeed so far to convert the short-chain expansion of the distribution 
function to a close expression. 

3. 3. The scattering function 

The scattering function of a semiflexible polymer is defined by 

S{q, ^) = ]^ ds2 dsr (exp(zq(r(s2) - v{s{)))) . (14) 

Expressing r(s2) — r(si) in (fT^ through the tangent vectors, r(s2) — r{si) = J^^ dst{s), 
and representing the average in ()14|) using the formalism of the quantum rigid rotator 
yields that the scattering function of the semiflexible polymer S{q, N) is the inverse 
Laplace transform of G{q,p)/p^ multiplied with the factor 2/N |21],[2ni- Thus, the 
short-chain expansion of G{k,p) enables one to derive in a straightforward way the 
short-chain expansion of the scattering function. Taking into account the first four 
terms in the sum over s in |TT|l results in the following expression for the scattering 
function 

,,,,,, X,,. 2A^ A^sinx lA^cosx 7 N'^smx 6 N'^ 

S4., N)/N = SAx)/N + -_-_ + -_ + _ —— + 

UN^sinx A N^cosx GN^cosx 32 N^ 307 A^3gin3, 



15 x3 15 x2 5 x^ 63 x^ 7560 x 
125iV3sinx A^^sinx 31 A^^cosx 157A^^cosx 31 

^ 1-3 = 1 t: N cosx, 

126 x3 x5 126 x2 63 x^ 7560 

where x = qN, and 

S,(x)//V = 2^-1+251M 

x^ x^ X 
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is the scattering function of a stiff rod. The plot of the scattering function multiphed 
by q for different values of Smax is shown in FigO The accuracy of the computations is 

3 



2 

z 

cr 

M 1- 
cr 

0- 

4 8 12 

qN 




Figure 1. The plot of qS{q,N) for chain length N = 1.5. Dots: stiff rod; continuous: 
Smax = 5; dashes: Smax = 8. 

determined by the values of x = qN, where curves corresponding to different Smax begin 
to diverge. The FigQ shows that for = 1.5 the continuous and dashed curves begin 
to diverge for x > 9. 

3.4- The extension-force relation 

The partition function of a semiflexible polymer with one end fixed and force f applied 
to the another end 



f 



N 



Z{f,N) = {exp{-^ dst{s))j (15) 

can be expressed through the distribution function of the end-to-end distance as follows 

Z{f,N)=Gik=-t-^,N). (16) 

Using the definition of the free energy F = —ksT In Z{f, N) the extension-force relation 
can be expressed through the partition function as 

dF _ d\nZif,N) 
^--Jf- df • (17) 

Thus, the short-chain expansion of the extension-force relation is directly obtained from 

that of the Fourier transform of the distribution function of the end-to-end distance 

fll2|) . Taking into account the first Smax = 4 terms in the series in (fTT|) results in the 

following expression for the extension-force relation 

R \ 11 1 1 

— = h coth X Nx coth x H — Nx coth^ x -\ A^^ coth x 

N X 6 6 6 90 

11 7 1 

H N'^x N'^x^ coth X N'^x coth^ x H N'^x^ coth^ x 

20 36 180 36 

1 11 19 1 

coth X H N^x^ coth x xN^ N^x^ coth^ x 



1512 1512 2520 120 
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N^x^ coth^ X + 



1 11 
N^x^ coth^ X + 



xcoth X, 



840 1890 216 2520 

where x = f'N. The first two terms on the right-hand side of ()18|) gives the extension- 
force relation for stiff rod. The log- log plot of 1 — R/N for different Smax and for stiff 
rod is shown in Fig|21 The slope for the stiff rod is —1, while the slope for the worm- like 



1,00 



0,37, 



0,14 



0,05 




0,14 0,37 1,00 2,72 
fN 



7,39 



Figure 2. Extension- force relation for chain length N = 1.1. Dots: stiff rod; 



continuous: s„ 



4; dashes: Sn 



6. 



chain approaches the value —1/2, which is the asymptotic result for finite worm-like 
chain at large forces [HO] . Fig|21 shows that the convergence of the short-chain expansion 
for the extension-force relation is worse than that of the scattering function. 

4. Conclusions 

To conclude, using the exact computation of a large number of moments of the end- 
to-end distribution function G{r, N) of the worm-like chain, which are obtained from 
the representation of the distribution function as a matrix element of the infinite order 
matrix [211,121]; we have established the analytical form of the coefficients in Taylor 
expansions of the moments for short A^. The knowledge of these coefficients enabled us 
to resummate the moment expansion of G{r, N) by taking into account consecutively 
the deviations of the moments from their stiff rod limit. Within this procedure we 
have derived the short-chain expansion for the distribution function of the end-to-end 
polymer distance, the scattering function, and the extension-force relation, by taking 
into account the deviations of the moments from their stiff rod limit to the seventh order 
in A^. The procedure can be extended to higher orders. The short-chain expansion could 
be useful in studies of the behaviour of short polymers, where the deviation from stiff 
rod is small. 
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